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ABSTRACT 

Aims. The goal of this work is to characterize the millimeter and centimeter properties of intermediate/high-mass young stellar objects 
(YSOs) in order to search for any evolutionary trend. 

Methods. We carried out observations at 1.2 mm with the IRAM 30 m telescope, and at 3.6 and 1.3 cm with the VLA toward a sample 
of 11 luminous (> 10 3 Lq) IRAS sources classified as high-mass protostellar object candidates. The most promising regions were 
additionally observed at 7 mm with the VLA. 

Results. The 1.2 mm emission is detected toward all the IRAS sources. In some cases the emission shows a clear peak surrounded by 
some substructure, while in other cases the emission is very extended and weak. The 1.2 mm emission is likely tracing the dust core 
in which the intermediate/high-mass YSO is forming, for which we estimated masses ranging from 10 to 140 Mq. For all (but one 
of) the sources, we detected centimeter emission associated with the IRAS source, being compact or ultra-compact (UC) H n region 
candidates. From the spectral index between 3.6 and 1.3 cm, we identified the possible emission mechanism in the centimeter range. 
Deconvolved sizes of the centimeter sources range from < 0.01 to 0.3 pc, and the physical parameters of the UCHn regions indicate 
that the ionizing stars are early B-type stars. The 7 mm emission is partially resolved for the four regions observed at this wavelength, 
and we estimated the contribution of the dust emission to the 7 mm flux density, ranging from negligible to 44%. By combining 
our data with the infrared surveys 2MASS, MSX and IRAS, and IRAC-Spitzer data (when available), we built the spectral energy 
distribution and fitted a modified blackbody law taking into account the contribution of the free-free centimeter emission. We found 
dust temperatures between 25 and 35 K, dust emissivity indices between 1.5 and 2.2, and masses similar to the masses derived from 
the 1.2 mm continuum emission. In addition, we found a correlation between the degree of disruption of the natal cloud, estimated 
from the fraction of dust emission associated with the centimeter source relative to the total amount of dust in its surroundings, and 
the size of the centimeter source. 

Conclusions. From the correlation found, we established an evolutionary sequence in which sources with compact millimeter emission 
clearly associated with compact centimeter emission are younger than sources with the millimeter emission dispersed and with the 
centimeter emission extended. Such a sequence is consistent with the evolutionary stage expected from maser/outflow/dense gas 
emission reported in the literature, and with the infrared excess of the 2MASS sources associated with the centimeter source, estimated 
in this work. 

Key words. Stars: formation - ISM: dust - ISM: H n regions - Radio continuum: ISM 

1 . Introduction i Patel et al.ll2005l). a nd infall motions (e. g.. lBeltran et al.ll2006ah 

[Zapata et alj I2008) have been detected in massive protostars, 
High-mass stars (M > 10 Mq) play a major role in the physi- sug g estin g that there is disk-mediated accretion in high-mass 
cal and chemical evolution of their host galaxies, through strong starSi and that they may form as a sca i e d-up version of the low- 
ultraviolet radiation and winds, and through supernova explo- mass stars HoweV er, the different evolutionary stages describing 
sions. However, the formation scenario of these stars is not how the massive protostar approaches the main-sequence may 
well understood, due to the large distances (> 1 kpc) and differ from those for i ow . mass stars> sinc e massive stars evolve 
heavy extinctions (A v > 100) of the regions where they form. much faster to the main . se quence (e. g., iBernasconi & Maided 
This, coupled with the fact that high-mass stars seem to form rf 996)i and? as a conse quence, the massive protostar produces 
in clustered mode (Testi et aH [1998] l2000t | Hillenbrand | [19951 TJV photons while is still accreting matter. Thus, the centimeter 
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[ Hillenbrand & Hartmann| | 1998[) , where a single massive dense emission, tracing the ionized gas, and the millimeter emission, 



core of gas forms a group of stars, makes observationally usua Uy tracing the dust, are crucial to understand the first stages 

more difficult to study these regions. A result recently found m tne f orma ti n of a massive star 
in massive s tar formation is that outflows, some of them colli- 
mated (e. g., iBeufher et al1l2002t iGibb et all 2003; Zh ang et 

ID The work presented here is part of a large project aimed 



12005b . disk-like structures (e. g., ICesaroni etail 1 19971 Il999t at studying in detail (with high sensitivity and high angu 

lar resolution) a sample of high-mass star-forming regions 

Send offprint requests to: Alvaro Sanchez-Monge, (ibol £ 1000 Lq) in different evolutionary stages, in order 

e-mail: asanchez@am.ub.es to characterize the millimeter and centimeter properties as the 
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Table 1. List of regions observed in this work 







Dist. 






Region 


(Lq) 


(kpc) 


Instrument 3 


Ket. 


IRAS 00117+6412 


1300 


1.8 


V + I 


1 


IRAS 04579+4703 


4000 


2.5 


V + I 


j 


IRAS 18123-1203 


7900 


3.1 


V + I 


1 


IRAS 18171-1548 b 


7900 


2.6 


V + I 




IRAS 18212-1320 


1600 


2.2 


I 




IRAS 19045+0813 


1600 


1.6 


V + I 




IRAS 20293+4007 


7900 


3.4 


V + I 




IRAS 22187+5559 


7900 


2.9 


V + I 




IRAS 22198+6336 


1300 


1.3 


V 




NGC 7538-IRS9 c 


40000 


2.8 


V 


2 


IRAS 23448+6010 


2500 


2.0 


V + I 


1 



a V: VLA, I: IRAM 30 m. 

This source appears as IRAS 18172-1548 in Molinari et al. 

c Identified wit h IRAS 23 118+61 10 

References: (D lMolinari et al.ll 19961 . (2) IMueller et~alll2002l . 



massive young s tellar object (YSO) e v olves. For this we se - 
lect ed, from the IMolinari et all d 19961) , iMueller et all (|2002|) . 
and lSridharan et alj d2002l) samples, regions located at distances 

< 3 kpc, in order to be sensitive down to masses < 1 Mq 
(from the millimeter emission), and to have a spatial resolution 

< 5000 AU. For each region, we searched the literature for cen- 
timeter and millimeter emission, and found 1 1 regions with no 
imaging at centimeter or millimeter observations. To complete 
the sample of nearby and luminous massive star-forming regions 
observed at millimeter and centimeter wavelengths, we carried 
out VLA and IRAM 30 m observations toward these 1 1 regions. 
Nine regions were observed in the millimeter range, and 10 in 
the centimeter range (8 coinciding with the regions observed at 
millimeter wavelengths, see Table[TJi. 

In § |2l we summarize the observations, in § [3] we give the 
general results for the millimeter and centimeter data, in §|4]we 
discuss the nature of the centimeter emission, comment the re- 
sults for each individual source, and discuss global properties 
possibly indicating different evolutionary stages of the massive 
YSOs. Finally, in §|5]we summarize our main conclusions. 

2. Observations 

The MPIfR 117-element bolometer array MAMBO-2 at the 
IRAM 30 m telescop43 was used to map the 1.2 mm dust contin- 
uum emission toward the 9 selected high-mass protostar candi- 
dates (see Table [TJ. The observations were carried out on 2004 
January 23 to 25. The main beam of the IRAM 30m at 1.2 mm 
has a HPBW size of 11". We used the on-the-fly mapping mode, 
in which the telescope scans in azimuth along each row in such 
a way that all pixels see the source once. The sampled area was 
typically 320" x 80", and the scanning speed was 6" s . Each 
scan was separated by 8" in elevation. The secondary mirror was 
wobbling at a rate of 2 Hz in azimuth with a wobbler throw of 
70". The average zenith opacity at 225 GHz was 0.3-0.4. The 
overall uncertainty in calibration is estimated to be 20%. Data 
reduction was performed with the MOPSIC software package 
(distributed by R. Zylka). 

1 IRAM is supported by INSU/CNRS (France), MPG (Germany), 
and IGN (Spain). 



The centimeter observations were carried out on 2004 July 

10 and 13 using the Very Large Array (VLA) of the NRAC0 in 
the D configuration. The regions observed at 3.6 cm, the phase 
calibrators used and their bootstrapped fluxes, as well as the pa- 
ramete rs of the synthesized beam, the robust parameter of Briggs 
(119951) . and the rms noise of the maps are given in Table [2] and 
the same information for 1.3 cm is given in Table [3] The posi- 
ti ons of the phas e cente rs c orrespond to the coo rdinates given 
in IMolinari etaH d 19961) or IMueller et ail d2002l) . In all cases, 
except for IRAS 19045+0813 and NGC 7538-IRS9, they corre- 
spond to the catalog positions of the IRAS sources. The inte- 
gration time for each source was about 5 minutes at 3.6 cm and 

11 minutes at 1.3 cm. Absolute flux calibration was achieved 
by observing 3C48, with an adopted flux density of 3.15 Jy and 
1.12 Jy at 3.6 and 1.3 cm, and 3C286, with 5.21 Jy and 2.52 Jy 
at 3.6 and 1.3 cm respectively. Images were constructed using 
a variety of weighting schemes by varying the Briggs robust pa- 
rameter from + 1 to +5 (roughly equivalent to natural weight) for 
most of the regions (see Table [2] and [3]). In order to increase the 
SNR of the faint sources at 1.3 cm, we tapered the Mv-data at 
50 kA, except for IRAS 22187+5559, for which we tapered the 
Mv-data at 25 kA, in order to recover the extended emission. 

The 7 mm observations carried out toward four regions (see 
Table |H were made with the VLA in the D configuration during 
the first months of 2007, with 8 EVLA antennas in the array. In 
order to minimize the effects of atmospheri c phase fluctuations, 
we u sed the technique of fast switching dCarilli & Holdawavl 
119971) between the source and the phase calibrator over a cy- 
cle of 120 seconds, with 80 seconds spent on the target and 40 
seconds on the calibrator. The integration time for each source 
was about 50 minutes, except for IRAS 04579+4703 with an on- 
source time of about 100 minutes. The images were constructed 
using natural weighting and tapering the Mv-data at 60 kA, except 
for IRAS 22198+6336, to increase the SNR. 

For simplicity, we will refer to each IRAS source by an "I" 
followed by the RA identification of the catalog name. 

3. Results 

3.1. Results at 1.2 mm 

In Table [5] we summarize the main results of the IRAM 30 m 
observations, in particular the peak intensity, flux density, and 
mass of gas and dust for each region. In the bottom panel of 
figures [JJ to [10] we show in grey scale and black contours the 
1.2 mm continuum emission. 

None of the regions observed with the IRAM 30 m shows 
a point-like morphology. On the contrary, four regions (118123, 
118171, 120293, and 123448) have very extended emission, with 
no clear peak. The other regions have a clear emission peak and 
show some substructure. For 118212 no millimeter source was 
detected (< 620 mJy). Thus, this source was not observed with 
the VLA. 

3.2. Results at 3.6 cm, 1.3 cm, and 7 mm 

In Tables [6] [7] and [8] we list all the sources detected in each 
region above the 5<x level at 3.6 cm, 1.3 cm, and 7 mm respec- 
tively. The three tables give, for each centimeter source identi- 
fied in each region (labeled in increasing RA), the coordinates, 



2 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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Table 2. Parameters of the observations at 3.6 cm with the VLA 
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Phase Center 


Phase 


Bootstrapped Flux 


Beam 


PA. 




rms 


R pen on 


a(J2000.0) 


<5(J2000.0) 


r^al lhratnr 

V Clll L'l Ll 1 1 '1 


(W) 


l ill V L 1 


\ ) 




( //Tv hpam ' ^ 

\Lijy L'^-iiiii i 


IRAS 00117+6412 


on \a 11 11 

UU l^f z / . / j 


_i_A/1 1Q A(\ 1 
TU i + Zo ^O.Z 


0102+584 


1.885 +0.003 


US x 8.5 


76 


+ 1 


38 


IRAS 04579+4703 


05 01 39.73 


+47 07 23.1 


0423+418 


1.986 ±0.004 


10.7 x 9.4 


48 


+5 


22 


IRAS 18123-1203 


18 15 07.33 


-12 02 42.0 


1822-096 


1.364 ±0.002 


12.3 x 8.1 


5 


+ 1 


43 


IRAS 18171-1548 


18 20 04.62 


-15 46 46.6 


1822-096 


1.364 ±0.002 


13.0x8.1 


2 


+ 1 


185 


IRAS 19045+0813 


19 06 59.77 


+08 18 42.8 


1922+155 


0.666 ± 0.002 


8.9x7.8 


43 


-3 


54 


IRAS 20293+4007 


20 31 07.92 


+40 17 22.8 


2025+337 


1.826 ±0.003 


7.5 x 6.5 


-37 


-5 


58 


IRAS 22187+5559 


22 20 34.89 


+56 14 39.5 


2148+611 


0.971 ±0.003 


7.8x6.6 


-1 


-5 


155 


IRAS 22198+6336 


22 21 27.62 


+63 51 42.2 


2148+611 


0.971 ±0.003 


10.4x8.8 


7 


+5 


27 


NGC 7538-IRS9 


23 14 01.68 


+61 27 19.9 


0014+612 


1.343 ±0.004 


10.2x9.3 


-14 


+5 


700 a 


IRAS 23448+6010 


23 47 20.07 


+60 27 21.1 


0102+584 


1.753 ±0.003 


10.0 x 8.7 


4 


+5 


30 



a rms limited by dynamic range. 



Table 3. Parameters of the observations at 1.3 cm with the VLA 





Phase Center 


Phase 


Bootstrapped Flux 


Beam 


PA. 




rms 


Region 


a(J2000.0) 


<5(J2000.0) 


Calibrator 


(Jy) 


(arcsec) 


C) 


Robust 


(/Jy beam 


IRAS 00117+6412 


00 14 27.73 


+64 28 46.2 


0102+584 


2.256 ±0.015 


5.9x4.5 


-21 


+5 


67 


IRAS 04579+4703 


05 01 39.73 


+47 07 23.1 


0423+418 


1.617 ±0.017 


5.4 x 4.4 


45 


+5 


85 


IRAS 18123-1203 


18 15 07.33 


-12 02 42.0 


1832-105 


0.913 ±0.015 


6.9x4.6 


-20 


+5 


63 


IRAS 18171-1548 


18 20 04.62 


-15 46 46.6 


1832-105 


0.913 ±0.015 


8.3 x 5.0 


-34 


+5 


112 


IRAS 19045+0813 


19 06 59.77 


+08 18 42.8 


1922+155 


0.504 ± 0.002 


6.0x4.2 


41 


+4 


65 


IRAS 20293+4007 


20 31 07.92 


+40 17 22.8 


2025+337 


1.946 ±0.013 


5.1 x 4.3 


41 


+5 


80 


IRAS 22187+5559 


22 20 34.89 


+56 14 39.5 


2137+510 


0.626 ± 0.008 


7.9x6.5 


50 


+5 


60 


IRAS 22198+6336 


22 21 27.62 


+63 51 42.2 


2230+697 


0.601 ± 0.005 


3.6x2.9 


-5 


+ 1 


63 


NGC 7538-IRS9 


23 14 01.68 


+61 27 19.9 


2230+697 


0.601 ± 0.005 


3.7 x 3.1 


14 


-5 


204 


IRAS 23448+010 


23 47 20.07 


+60 27 21.1 


0102+584 


2.256 ±0.015 


5.3x4.3 


40 


+5 


76 



Table 4. Parameters of the observations at 7 mm with the VLA 





Phase Center 


Phase 


Bootstrapped Flux 


Beam 


PA. 




rms 


Region 


a(J2000.0) 


5(J2000.0) 


Calibrator 


(Jy) 


(arcsec) 


C) 


Robust 


(/Jy bearrT 1 ) 


IRAS 04579+4703 


05 01 39.73 


+47 07 23.1 


0502+416 


0.424 ±0.012 


3.7x3.2 


40 


+5 


94 


IRAS 19045+0813 


19 06 59.77 


+08 18 42.8 


1856+061 


0.249 ± 0.006 


3.9x3.3 


44 


+5 


100 


IRAS 22187+5559 


22 20 34.89 


+56 14 39.5 


2250+558 


0.300 ± 0.019 


3.6x3.1 


49 


+5 


178 


IRAS 22198+6336 


22 21 27.62 


+63 51 42.2 


2148+611 


0.530 ±0.032 


2.2 x 1.3 


-43 


+0 


226 



peak intensity, flux density, and deconvolved size and position 
angle obtained by fitting a gaussian with the task JMFIT in AIPS. 
In all three tables, we classify a VLA source as "central" if it 
is located within the central 50" of diameter around the IRAS 
source, corresponding to 0.5 pc at an average distance of 2 kpc, 
this being the expected size of the cluster forming around a high - 
mass source (iHUlenbrand & Hartmanr][l99l iTesti et al.| [l999). 
For sources with extended emission, we estimated the intensity 
peak and flux density using the TVSTAT task in AIPS. 

In the case of NGC 7538-IRS9 at 3.6 cm, the rms was limited 
by dynamic range from the strong source VLA 1 (NGC 7538- 
IRS1) and the extended emission toward the north of VLA 1 
(which corresponds to an Hn region detected at 21 cm, see 
Fig. [9]). We performed self-calibration in phase by using the 
strong source VLA 1. Subsequently, we substracted the clean 
components of VLA 1 using the task UVSUB in AIPS and used 
different wv-ranges in order to better clean the extended emission 
and reduce the rms noise. Although the rms noise was lowered 



from 1.28 to 0.7 mJy beam 1 , this was not low en ough to detect 
VLA 2 and VLA 3, for which lSandell et all d2005l) report a flux 
density of 0.76 + 0.15 mJy. 

We estimated the spectral index a (S v oc v") of the central 
sources, listed in Table [9] using the flux densities at 3.6 and 
1.3 cm. In order to properly compare the flux at both wave- 
lengths, we made 3.6 and 1 .3 cm images using the same wv-range 
(4-25 k/l) at both wavelengths. For the sources not detected at 
both wavelengths in the 4-25 k/l maps, we used the flux densities 
obtained with no restrictions in the wv-range, listed in Tables|6]or 
Q The upper-limit for non-detections was set to 4 times the rms 
noise level of the map, or slightly higher for extended sources. 
For details in the calculation of the spectral indices see the foot- 
notes of Table [9] 

Addi tionally, we search ed the NRAO VLA Sky Survey 
(NVSS, ICondon et afl [l998) for emission at 21 cm, and esti- 
mated the spectral index between 21 and 3.6 cm. We made new 
images at 3.6 cm with an Mv-range of projected baselines of 1- 
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Table 5. Parameters of the 
IRAM 30 m 


1 .2 mm sources 


observed with 




rms 


r peak 

iy 




Mass b 


Region 


(mJy beam - 


') (mJybeam -1 ) 


(mJy) 


(Mq) 


IRAS 00117+6412 


8 


108 


1200 


36 


IRAS 04579+4703 


7 


92 


390 


23 


IRAS 18123-1203 


15 c 


68 


440 


41 


IRAS 18171-1548 


16 c 


74 


2200 


142 


IRAS 18212-1320 


14 c 


<56 


< 620 




IRAS 19045+0813 


9 


100 


390 


10 


IRAS 20293+4007 


9 


55 


570 


64 


IRAS 22187+5559 










VLA 3 d 


7 


46 


180 


14 


VLA5 d 


7 


86 


390 


32 


IRAS 23448+6010 


11 


63 


1600 


63 



a Flux density inside an aperture of 80" of diameter, except for IRAS 

22187+5559, where the aperture was 50" for the source associated with 

VLA 5 and 30" for the source associated with VLA 3. 
b Mass of gas and dust derived assumi ng T& = 30 K, and a dust mas s opacity 

coefficient of 0.9 cm 2 g _1 at 1.2 mm lOssenkopf & Henning 1994). 
High rms noise due to the low elevation of the sources at the time of the 

observations. 
d See Section l4~2~7l 



5 kA, which corresponds to the wv-range of the NVSS images. In 
case of non-detection we followed the same procedure described 
above. 

For the four regions additionally observed at 7 mm (104579, 
119045, 122187, and 122198), we detected emission in all the 
cases (see Table[8]l, with 119045 being the faintest source. For all 
the sources, the 7 mm emission was partially resolved. We esti- 
mated the contribution of ionized gas to the 7 mm emission by 
using the spectral index from 3.6 to 1.3 cm (Tabled and extrap- 
olating the flux density at 1 .3 cm (Tabled to 7 mm. From the es- 
timated free-free flux density at 7 mm and the total observed flux 
density (from Table |H), we estimated the contribution from dust 
at 7 mm, which was - 38% for the case of 104579, ~ 44% for 
122198, and negligible for 119045 and 122187 (Table[TT)- We ad- 
ditionally estimated the mass of the dust compone nt at 7 mm as- 
sumin g a dust temperature in the 50-100 K range (Gome zet alj 
120031) . a dust emissivity index of 2 and a dust mass opacity coef - 
ficient of 0.9 cm 2 g _1 at 1.2 mm (lOssenkopf & Henningll 1994b . 
Note that we assumed a higher temperature for the dust emission 
at 7 mm than for the IRAM 30 m emission at 1 .2 mm, because 
the IRAM 30 m emission must come mainly from the (large- 
scale) envelope of the massive object, while the 7 mm emis- 
sion is tracing the compact dust component possibly related with 
a disk, which must be hotter than the envelope. The resulting 
masses for the 7 mm dust components are listed in Table [TT] 

Figures [TJ to [10] show the maps for each region. For most of 
the regions there are two main panels showing the /T-band of 
the Two Micron All Sky Survey (2MASS), and the IRAM 30 m 

1.2 mm continuum emission, together with the VLA 3.6 and 

1.3 cm continuum emission, and 7 mm continuum emission 
when available. For regions with sources beyond the central 
~ 50" around the IRAS source (Figs. [3][5j and|6]), there is an 
extra top panel showing the area covered by the primary beam at 
3.6 cm (~5'), while for Figs. [7] [8] and|9] the distribution of the 
panels is slightly different. 



3.3. Infrared photometry from 2MASS and /RAC-Spitzer 

For each region, we downloaded the sources of the 2MASS 
catalog within the VLA primary beam at 3.6 cm, and we es- 
timated for each 2MASS source the infrared excess from the 
(/ - H) vs (H - K) diagram, measured as the difference be- 
tween the (H - K) color and the (H - K) color corresponding 
to a reddened main-sequen ce star (following the reddening law 
of lRieke & L ebof sky 1985b . We assigned a color and a different 
star symbol to each infrared excess interval: three-points blue 
stars for infrared excesses < -0.4; four-points green stars for 
the interval (-0.4-0), corresponding to main-sequence stars, gi- 
ants, and Class III and II sources with small infrared excess; 
five-points yellow stars for (0-0.4), corresponding to Class II 
sources; six-points orange stars for (0.4-1), corresponding to 
Class I sources, and ten-poin ts red stars for infrared excesses 
> 1 (Matsuyanagi et al. 2006). The 2MASS sources are marked 
as stars with the corresponding color in Fig.[TJto Fig.QJJ In our 
sample, 104579 and NGC 7538-IRS9 are the two regions where 
the 2MASS source associated with the IRAS source shows the 
largest infrared excess. Large infrared excesses are interpreted 
as reflecting an early evolutionary stage not only for low- mass 
YSOs (e. g.. lStrom et al.ll 19931: [Renvoi] i & Hartmann II 19951). but 
also for interme diate/high-mass YSOs (e. g.. iMassi et alJll9"99l 
lNurnbergerl2"003l) . 

In addition, 118171 and 119045 were observed with IR AC- 
Spitzer, whose images are shown in lKumar & Gravel (2007). For 
these two regions, we studied the IRAC sources within the cen- 
tral 120" around the UCHii region, and plotted the sources in 
#-[3.6] vs [3.6]-[4.5] and [3.6]-[4.5] vs [4.5]-[5.8] diagrams, 
allowing us to distinguish between Class II/III and Class 0/1 
sourc es d Allen et al J 120041: iHartmann etail 120051: iGetman et all 
120071) . For simplicity, in Figs. H to we plotted only the IRAC 
sources classified as Class 0/1 sources and with no 2MASS coun- 
terpart. We contrasted the classification of YSOs from IRAC col- 
ors with the classification from JHK infrared excess (described 
above) and found that in general IRAC and 2MASS colors yield 
a similar classification. Thus, for the regions not observed with 
IRAC we consider the classification of YSOs from JHK infrared 
excess as a first valid approach to the evolutionary stage of each 
source. 

3.4. Spectral Energy Distributions 

We searched the liter ature for data from 2MASS, IRAC-Spitzer, 
MSX, and SCUBA dDi Francesco et alj|2008l) for each region. 
We found an infrared counterpart of the IRAS source for all the 
cases except for 122198. With these gathered data, and the data 
presented in this work, we built the SED for each IRAS source, 
and fitted a modified blackbody law for the far-infrared to mil- 
limeter data, taking into account the contribution of a free-free 
optically thin law for the centimeter data. Note that we did not 
attempt to fit the infrared excess coming from components with 
higher temperatures. We ado pted a dust mass opacity coefficient 
at 1.2 mm of 0.9 cm g - 1 (lOssenkopf & Henning|[T994T) . a gas- 
to-dust mass ratio of 100, and assumed that the opacity follows 
a power law. The modified blackbody law allowed us to esti- 
mate a range of masses, dust temperatures, and dust emissivity 
indices associated with the envelope of each intermediate/high- 
mass YSO (see Table [TTI and Fig. [TT}. The median values were 
~ 25 Mq, ~ 29 K, and 1.9, respectively. Dust temperatures and 
dust emissivity indices obtain e d in th is work are similar to those 
obtained by iMookeriea et alj d2007l) toward a sample of high- 
mass protostellar objects located in the outer Galaxy. For 120293 
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Table 6. Parameters of the sources detected at 3.6 cm 



5 





VLA 


Central 


o-(J2000.0) 


<5(J2000.0) 


.peak h 

1 V 


b 

Sy 


Deconv.Size 


RA. 


Region 


Source 


Source 11 


( hm 5 ) 


( ° ' " ) 


(mJy beam" 1 ) 


(mJy) 


(") 


O 


IRAS 00117+6412 


2-3 c 


Y 


00 14 28.12 


+64 28 46.2 


1.23 ±0.04 


1.61 ±0.08 


8.8 x 4.3 


73 


TD AC Ozl^'7Gj-zl7fn 


i 
i 


v 


ni en 

Uj ui jy.yz 


_i_ An fi7 9 1 1 
++ / u / Zl . 1 


U. JO ± U.Uj 


fiKj.fi C\A 
U. 1 J ± U.LH- 


J.J x u.u 


jU 


TP A ^ 1 S1 1^ 1 9fH 


i 
i 


IN 


is 1 s m ^ 

1(5 1J UI.jO 


19 01 zL^ 9 
— 1Z Ul HO.Z 


u.z / ± u.uo 


fiCC j.fl7'? 
U.OO ± U.ZJ 


1 Q 9 v 1 1 1 

A 11.1 


1 M 

IO J 




2 


Y 


18 15 07.55 


-12 01 43.3 


0.22 ± 0.04 


0.45 ±0.12 


12.2 x 8.1 


19 




3 


N 


18 15 08.01 


-12 05 10.0 


2.25 ± 0.09 d 


2.62 ± 0.17 d 


5.4 x 2.9 


2 




4 


Y 


18 15 08.05 


-12 01 58.0 


0.23 + 0.04 


0.53 ±0.14 


12.4 x 9.5 


107 




5 


Y 


18 15 08.12 


-12 01 38.0 


0.30 + 0.04 


0.58 ±0.12 


10.4 x 8.2 


81 




6 


Y 


18 15 08.75 


-12 01 58.0 


0.23 ± 0.04 


0.43 ±0.12 


10.9 x 6.1 


104 




7 


N 


18 15 14.42 


-11 59 32.0 


0.79 ± 0.16 d 


1.85 ± 0.55 d 


15.2 x 8.0 


52 


IRAS 1 8 17 1 — 1548 


1 


Y 


1 8 20 06.40 


1 r- a S* A A 1 

-15 46 44.1 


1.57 ± 0.21 


3.42 ± 0.63 


12.2 x 9.9 


18 


IRAS 19045+0813 


1 


N 


19 06 50.61 


+08 19 20.6 


0.50 ± 0.05 d 


0.54 ± 0.09 d 


4.1 x 0.0 


64 




2 


Y 


19 06 58.64 


+08 18 58.3 


0.23 ± 0.05 


0.38 ±0.11 


16.2 x 0.0 


21 




3 


Y 


19 06 59.34 


+08 18 56.7 


0.41 ± 0.04 


0.97 ± 0.14 


12.3 x 7.2 


91 




4 


Y 


19 06 59.63 


+08 19 08.9 


0.54 + 0.04 


1.41 ± 0.15 


15.5 x 6.4 


48 




J 


v 
i 


1 Q HA in 
IV Uo jy. 1 1 


, no in in I 
+Uo ly Ly.l 


V.jV + U.U4 


U. JZ ± u.l 1 


in 1 v 9/i 

1U.J X Z.4 


190 




6 


Y 


19 07 00.12 


+08 19 09.9 


0.50 ± 0.04 


1.32 ± 0.15 


13.2 x 7.7 


130 




2-6 c 


Y 






0.56 + 0.13 


3.14 ± 0.19 








7 


Y 


19 07 00.45 


+08 18 44.5 


0.59 + 0.05 


0.60 ± 0.08 


3.4 x 0.0 


152 


IRAS 20293+4007 


1 


Y 


20 31 07.17 


+40 17 20.9 


0.87 ± 0.06 


2.84 ± 0.24 


12.4 x 8.7 


22 




2 


Y 


20 31 07.19 


a r\ 1 n 1 

+40 17 27.1 


0.72 ± 0.06 


2.08 ± 0.22 


10.6 x 8.7 


179 




1 9C 
1— Z 


v 
i 






n on _i_ n 11 


■3 7Q j_ fl Q 1 

j. ly ± U. ji 








3 


N 


20 31 19.11 


+40 18 09.7 


1.39 ± 0.10 d 


1.63 ± 0.20 d 


4.4 x 0.0 


72 


TT1 AC Tl 1 O "7 i CCCn 

IRAS 2218/+5559 


1 


Y 


22 20 32.01 


+56 14 53.5 


1 cn 1 f\ n a 

1.59 ± 0.24 


1 "7 1 i A O A 

1. 11 ± 0.24 








2 


Y 


22 20 33.45 


+56 15 01.5 


2.24 ± 0.37 


3.31 ±0.37 








3 


Y 


22 20 33.67 


+56 14 26.4 


3.25 + 0.10 


32.58 ± 1.10 


39.0 x 10.7 


128 




4 


Y 


22 20 34.51 


+56 14 57.6 


3.21+0.10 


20.10 ±0.72 


19.3 x 13.9 


155 




5 


Y 


22 20 35.58 


+56 14 45.3 


5.59 ±0.10 


38.18 ±0.78 


21.3 x 14.0 


170 




1-5 C 


Y 






5.96 ± 1.06 


109 ±2 






IRAS 22198+6336 


1 


N 


22 21 24.48 


+63 52 20.2 


0.30 ± 0.02 


0.30 ± 0.03 


2.2 x 0.0 


26 




2 


Y 


22 21 26.68 


+63 51 38.2 


0.59 + 0.02 


0.57 ± 0.03 


1.5 x 0.0 


171 


NGC7538-IRS9 C 


1 


N 


23 13 45.47 


+61 28 19.7 


800 ± 10 d 


1360 ± 30 d 


10.2 x 6.0 


4 


IRAS 23448+6010 


1 


Y 


23 47 16.29 


+60 27 13.6 


0.39 ± 0.03 


0.56 ± 0.06 


13.3 x 0.0 


176 




2 


Y 


23 47 20.08 


+60 27 22.8 


0.40 + 0.03 


1.91 ±0.14 


23.7 x 13.4 


151 




4 


Y 


23 47 21.97 


+60 27 48.3 


0.21 + 0.03 


0.53 ± 0.09 


14.8 x 8.6 


78 




5 


Y 


23 47 22.47 


+60 27 12.7 


0.29 ± 0.03 


0.90 ±0.10 


15.0 x 11.7 


112 




6 


Y 


23 47 23.21 


+60 27 33.8 


0.26 + 0.03 


0.60 ± 0.08 


13.4 x 8.6 


7 




2-6 c 


Y 






0.40 ±0.10 


3.19 ±0.13 







11 Y = source within the central 50 " around the IRAS source; N = source out of the central 50 " around the IRAS source. 
b Intensity and flux density corrected for primary beam response. 
c Emission from the sources included in the range. 
d Primary beam correction factor greater than 1.5. 

c Extended emission is found toward the north of VLA 1, with an intensity peak of ~ 50 mJy beam -1 and a flux density of ~ 1.8 Jy. The 
coordinates of the intensity peak are a(J2000.0) = 23 h 13 m 32 s .63, and <5(J2000.0) = +61°29'53".7 (see Fig.|9]l. 



and 123448 we did not build the SED because the millimeter 
peak and the UCH n region (associated with the IRAS source for 
123448) are separated ~ 25" and 40", respectively, and could be 
tracing different objects. 

Regarding the contribution of the free-free emission, we 
adopted the spectral index listed in Table |9] (for NGC7538- 
IRS9 the spectral index used in the fit corresponds to -0.1). For 
122198 we downloaded 21 cm and 6 cm VLA continuum data 
from the archive (project AO 145) and reduced it. At the position 
of VLA 2 we found a compact source at 2 1 cm and 6 cm, with 
a flux density of 0.23 ± 0.05 mJy and 0.38 + 0.03 mJy, respec- 
tively. These values are in agreement with those obtained by us 
at 3.6 cm, 1.3 cm and 7 mm (see the SED in Fig.fTTl). 



4. Discussion 

4.1. On the nature of the cm emission associated with IRAS 
sources 

Most of the sources classified as central sources in Tables [6] 
and [7] have a flat or positive spectral index (see Table [9). All 
of them (with the exception of 118123) are nearly coincident 
with the corresponding IRAS source, which in all cases has a 
bolometric luminosity > 10 3 Lq, suggesting that the centime- 
ter sources associated with the IRAS source are tracing UCH n 
regions. However, the centimeter emission from the sources as- 
sociated with molecular outflows, 122198 and NGC 7538-IRS9, 
could also arise in gas ionized by shocks in the outflow. The cen- 
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Table 7. Parameters of the sources detected at 1.3 cm 
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( ° "' ) 


(mjy beam -1 ) 
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IRAS 00117+6412 


1 


N 


00 14 22.26 


+64 28 28.7 


0.47 ± 


0.08 


0.39 ±0.14 


unresolved 


- 




2 


Y 


00 14 28.32 


+64 28 41.2 


0.47 ± 


0.06 


0.70 ±0.15 


6.0 x 


0.0 


57 




3 


Y 


00 14 28.43 


+64 28 47.7 


0.57 ± 


0.06 


1.24 ± 0.19 


5.9 x 


5.2 


119 




2-3 c 


Y 






0.59 ± 


0.15 


1.72 ± 0.20 








IRAS 04579+4703 


1 


Y 


05 01 40.05 


+47 07 21.5 


0.54 ± 


0.09 


0.44 ±0.13 


2.2 x 


0.0 


31 


TRAS 1 8171 -1 548 


1 


Y 


1 8 20 06 46 


-1 5 46 39 9 


0.66 ± 


0.12 


2 77 + 62 


14.3 x 
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Tn A O 1 f\C\ A C . AO 1 O 

IRAS 19045+0813 


7 


Y 


19 07 00.59 


+08 18 44.1 


0.75 ± 


0.05 


1.10 ± 0.12 


5.3 x 


0.0 


143 


TT"> A C OAOAI i A AA^ 

1KAS 20293+400/ 


2 


Y 


TA 1 1 A"7 1 C 

20 31 0/.15 


i A A 1 1 1 1 

+40 1 / 28.1 


0.55 ± 


0.08 


1 AO i A Tl 

1.08 ± 0.22 


6.1 x 


3.1 


76 


IRAS 22187+5559 


1 


Y 


22 20 31.87 


+56 14 55.3 


0.70 ± 


0.05 


1.32 ±0.14 


8.4 x 


5.0 


12 




2 


Y 


22 20 33.38 


+56 15 06.8 


1.13 ± 


0.05 


3.79 ± 0.22 


13.2 x 


8.8 


179 






Y 


77 70 11 77 


+ Sn 14. 7Q 1 


1.99 ± 


0.05 


If.JJ X U.J1 


38.1 x 


8.3 


127 




4 


Y 


22 20 34.59 


+56 15 00.4 


1.95 + 


0.05 


6.05 ± 0.21 


11.2 x 


9.3 


99 




5 


Y 


22 20 35.64 


+56 14 46.4 


4.19 ± 


0.05 


19.31 ±0.28 


16.0 x 11.2 


170 




1-5 C 


Y 






4.64 ± 


0.78 


39.32 ± 1.34 








IRAS 22198+6336 


2 


Y 


22 21 26.87 


+63 51 37.2 


0.49 ± 


0.04 


1.11 ±0.13 


6.2 x 


1.4 


143 


NGC 7538-IRS9 


2 


Y 


23 14 01.10 


+61 27 18.8 


1.08 ± 


0.20 


1.43 ± 0.43 


4.6 x 


0.0 


13 




3 


Y 


23 14 01.79 


+61 27 20.0 


2.12 ± 


0.20 


3.83 ± 0.52 


4.4 x 


1.8 


14 


IRAS 23448+6010 


3 


Y 


23 47 20.61 


+60 27 16.0 


0.63 ± 


0.08 


0.80 ±0.15 


4.2 x 


0.0 


154 



11 Y = source within the central 50 " around the IRAS source; N = source out of the central 50 " around the IRAS source. 
b Intensity and flux density corrected for primary beam response. 
c Emission from the sources included in the range. 



Table 8. Parameters of the sources detected at 7 mm 
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Central 


or(J2000.0) 
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Deconv.Size 
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Source a 


("') 


( ° "' ) 


(mJy beam -1 ) 
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o 


IRAS 04579+4703 


1 


Y 


05 01 40.01 


+47 07 20.5 


0.67 ± 0.08 


1.45 ±0.26 


5.0 x 


2.5 


51 


IRAS 19045+0813 


7 


Y 


19 07 00.51 


+08 18 44.3 


0.58 ±0.11 


0.70 ± 0.20 


3.1 x 


0.0 


31 


IRAS 22187+5559 c 


3 


Y 


22 20 33.55 


+56 14 31.7 


1.09 ±0.28 


7.87 ±0.48 










5a 


Y 


22 20 35.27 


+56 14 36.2 


1.10 ±0.32 


2.81 ±0.51 


8.1 x 


2.7 


30 




5b 


Y 


22 20 35.56 


+56 14 45.4 


1.42 ±0.33 


3.02 ± 0.63 


6.6 x 


2.8 


110 




5c 


Y 


22 20 36.10 


+56 14 41.7 


0.99 ±0.16 


1.70 ±0.66 


6.7 x 


4.8 


5 




5d 


Y 


22 20 36.27 


+56 14 38.1 


1.07 ±0.25 


1.71 ±0.59 


4.6 x 


3.2 


140 




5a-5d d 


Y 






1.42 ±0.31 


12.05 ±0.48 








IRAS 22198+6336 


2 


Y 


22 21 26.76 


+63 51 37.9 


2.10 ±0.26 


2.76 ±0.53 


1.5 x 


0.3 


110 



11 Y = source within the central 50 " around the IRAS source. 

b Intensity and flux density corrected for primary beam response. 

c VLA 5a to VLA 5d are the subclumps found toward the position of VLA 5. 

d Emission from the sources included in the range. 



timeter continuum luminosity, S v d z , from shock-ionized gas is 
proportional to the outflow momentum rate, P, with a constant 
of proportionality including an efficiency fac tor, rj, defined as th e 
fraction of the stellar wind that is shocked (ICuriel et alJ fl987). 
For low-mass YSOs. I Angladal (119951) found that the centimeter 
emission can be produced by shock-ionized gas with efficiencies 
of 10%. In our case, the efficiencies needed to account for the 
observed centimeter emission should be much higher than those 
observed for low-mass YSOs, ~ 90% for 122198 and ~ 50% 
for NGC 7538-IRS9, being reasonable to assume that, at least 
for these two sources, the contribution of shock-ionized gas in 
these sources if any, is only a small fraction of the observed con- 
tinuum centimeter emission. Since these two sources seem to be 



among the youngest in our sample (see § 14. 3t , the contribution 
of shock-ionized gas in the other sources is expected to be also 
a small fraction of the observed centimeter emission. Thus, the 
centimeter sources associated with the IRAS source are most 
likely tracing compact or UCH n regions, and in Table [10] we 
list their physical parameters. Note that the emission measure, 
electron density, and opacity are beam-averaged values, and due 
to our angular resolution of ~ 10" at 3.6 cm, these values are 
affected by an important filling factor. 

From Table [TUl we can see that the underlying star exciting 
the UCH ii region is in all cases an early-type B star. We note 
that in the estimation of the physical parameters we assumed 
that the emission is optically thin. This is supported by the spec- 
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Table 9. Spectral indices of the central sources 





VLA 


Spectral Index" 


Spectral Index 




Region 


Source 


(zl— j.o cm) 


(i.o-i.j cm) 


Nature 


IRAS 00117+6412 


2-3 d 


> +0.13 


+0.18 + 0.20 


UCHu 


IRAS 04579+4703 


1 


> -0.2 


+ 1.1 ±0.4 C 


UCH ii / radiojet 


IRAS 18123-1203 


2 




< +0.9 C 






4 




< +0.8 e 






5 




< -0.4 C 






6 




< -0.T 






2-6 d 


> +0.17 






TR A S 18171 — 1 548 


1 


-~ KJ.J 




pnmnnpf 1-1 tt 


td a c i on/i ^ i no 1 1 


Z 




< +1.1 






3 




< +U.U 






4 




< -0.3 C 






5 




< +0.3 e 






6 




< -0.2 e 






2-6 d 


> -0.03 f 








7 


> +0.3 


+0.4 ±0.2 


UCH ii / radiojet 




i 




< — U.o 






2 




+0.1 +0.2 


CC\XX\V\5\C\ W TT 
l^UlllUuA^L 1111 




l-2 d 


> -0.2 






IRAS 22187+5559 


1 




-0.3 ±0.2 e 


ext. ionized? 




2 




+0.14±0.13 c 


ext. ionized? 




3 




+0.0 ±0.3 


compact H n 




4 




+0.6 ±0.5 


ext. ionized? 




5 




-0.06 ±0.17 


compact H n 




l-5 d 


+0.04 ± 0.02 


+0.13 ±0.09 




IRAS 22198+6336 


2 


+0.5 ±0.2 E 


+0.5 ±0.3 


UCH ii / radiojet 




Z 




> — l.Z 


UCH ii / radiojet 




3 




. f\ i e 

> +U.1 


UCH ii / radiojet 




2-3 d 


> -0.9 f 






IRAS 23448+6010 


1 




< -0.2 


UCH ii 




2 




< -0.2 


UCH ii 




3 




> -0.1 


UCHn 




4 




< +1.2 C 






5 




< +0.8 C 






6 




< +0.9 e 






l-6 d 


> +0.2 







11 Spectral index estimated from the flux density values at 21 cm from the NVSS iCondon et alJI 19981) and at 3.6 cm (present work) obtained 
with the i/v-range common to both wavelengths, 1-5 k/i. In case of non-detection at one of the frequencies, an upper limit of 4 times the rms 
noise of the map was used. 

b Spectral index estimated from the flux density values at 3.6 cm and 1.3 cm (present work) obtained with the wv-range common to both 

wavelengths, 4-25 kA. In case of non-detection at one of the frequencies, an upper limit of 4 times the rms noise of the map was used. 
c See main text. 

d Emission from the sources included in the range. 

e Non-detection at 3.6 cm and 1.3 cm in the same wv-range maps. Flux density at 3.6 cm or 1.3 cm from Tables[6]or|7] Upper limits at 3.6 cm 
or 1.3 cm (for the case of extended sources) as 4A 80 times the rms of the map (from Tables[2]or [3), with A = 1 + Q. s /£l h , where fl s is the solid 
angle of the source at the frequency where it is detected, and fib is the beam solid angle (see Beltran e t alj|200lb . 

f Non-detection at 21 cm and 3.6 cm in the same w-range m aps. Flux density at 3.6 cm for 119045 from Table[6] and flux density at 3.6 cm of 
0.76 ± 0. 15 mJv for NGC 7538-9 from lSandell et all d2005l) . 

° Flux density at 21 cm of 0.23 ± 0.05 mjy from the VLA archive data. See Section [3~4l 



tral indices close to zero reported in Table [9] for almost all the 
sources, except for 104579, for which the estimation of th e rate 
of ionizing photons is a lower limit (e. g.. lKeto et al.l l2008). and 
the spectral type of the ionizing star is earlier than B2-B3. For 
almost all the regions, the bolometric luminosity derived from 
the centimete r emission (from the spectral type and following 
Panagia 1973) is lower than the IRAS bolometric luminosity es- 
timated from the infrared emission. Only for 1001 17 and 122187 
we found the IRAS luminosity to be a factor of < 2 smaller. 



However, in both cases the centimeter sources split up in differ- 
ent components when observed with higher angular resolution 
(see § 14.2. 11 and |4.2.7t , being the 3.6 cm flux and the derived 
luminosity (at least for these two sources) an upper limit to the 
luminosity coming from the massive YSO. Thus, the IRAS lumi- 
nosity seems to be enough to account for the centimeter emission 
by photoionization. We concluded that all the centimeter sources 
associated with the IRAS sources of our sample seem to trace a 
compact or UCH n region. 
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Table 10. Physical parameters of the H n regions 





VLA 


V 


Diameter 




/£M a 








% 


Spectral 


Region 


Source 


(GHz) 


(pc) 


(K) 


(10 3 cm~ 6 pc) 


(10 3 cirr 3 ) 


(Mq) 




/ — 1 \ 

(s ') 


Type c 


IRAS 00117+6412 


2-3 


8.5 


0.054 


0.82 


22.3 


0.64 


1.3 x 10" 


3 


3.8 x 10 44 


B2 


IRAS 04579+4703 


1 


8.5 


< 0.040 


> 0.27 


> 7.2 


> 0.42 


< 1.9 x 10 


3 


6.8 x 10 43 


B2-B3 


IRAS 18171-1845 


1 


8.5 


0.139 


0.56 


14.8 


0.33 


1.8 x 10" 


2 


1.7 x 10 45 


Bl 


IRAS 19045+0813 


7 


8.5 


0.026 


1.02 


27.1 


1.01 


9.8 x 10" 


4 


1.1 x 10 44 


B2-B3 


IRAS 20293+4007 


2 


8.5 


0.158 


0.43 


11.8 


0.27 


1.9 x 10" 


2 


1.7 x 10 45 


Bl 


IRAS 22187+5559 


3 


8.5 


0.287 


1.53 


40.8 


0.38 


1.2 x 10" 




2.0 x 10 46 


B0.5 




5 


8.5 


n o a i 
U.Z4i 


2.51 


66.y 


U.53 


1.1 X 10 


1 


Z.j X 10 


BO. J 


IRAS 22198+6336 


2 


8.5 


< 0.010 


> 4.97 


> 132.0 


> 3.74 


< 6.4 x 10 


-4 


7.0 x 10 43 


B3 


NGC7538-IRS9 


2 


22.5 


< 0.062 


> 0.18 


> 39.1 


>0.79 


< 1.9 x 10 


-3 


9.0 x 10 44 


B1-B2 




3 


22.5 


0.038 


1.27 


279.0 


2.70 


4.0 x 10" 


-3 


2.4 x 10 45 


Bl 


IRAS 23448+6010 


1 


8.5 


< 0.129 


> 0.06 


> 1.7 


> 0.11 


< 3.1 x 10 


3 


1.6 x 10 44 


B2-B3 




2 


8.5 


0.173 


0.12 


3.1 


0.14 


9.0 x 10" 


-3 


5.5 x 10 44 


B2 




3 


22.5 


< 0.041 


> 0.12 


>25.9 


>0.80 


< 7.0 x 10 


-3 


2.6 x 10 44 


B2 



11 / corresponds to the filling factor. 

b Mass of ionize d gas estimated from the beam averaged electron density and the observed size of the source. 
c Estimated from Panagial dl973h . 



4.2. Comments on individual sources 
4.2.1. IRAS 00117+6412 

At 3.6 cm there is one source elongated roughly in the east- west 
direction, and at 1.3 cm the source splits up into two compo- 
nents, with the northern one, VLA 3, nearly coincident with the 
peak at 3.6 cm (Fig. |TJ. The millimeter emission has the main 
peak located ~ 15" to the west of VLA 3 and shows a rather cir- 
cular structure with some subcondensations. There are two H2O 
masers close to the m ain and the secondary millimeter peaks 
dCesaroni et al.|[l988h . suggesting that these two peaks are trac- 
ing embedded YSOs, and there is also out flow emission cen- 
tered around the sec ondary millimeter peak (Zhang et alj|2005t 
lKim&Kurd l2006). Observations with higher angular resolu- 
tion at 3.6 cm (Busquet et al., in prep.) reveal no emission at the 
position of VLA 2, and that the source detected in this work at 
3.6 cm consists of two components. One associated with VLA 3 
and with the brightest /T-band infrared source in the region, 
which has some infrared excess, and the other coincident with 
the main millimeter peak (likely harboring a deeply embedded 
intermediate-mass YSO). Thus, VLA 3 is most likely an UCH n 
region at the border of a dusty cloud where active star formation 
is taking place. The SED of 1001 17 can be fitted for temperatures 
between 26 and 30 K, and for masses between 21 and 30 Mq, 
which are among the median values of all the sources in this 
survey. 



4.2.2. IRAS 04579+4703 

The 3.6 cm source in this region is elongated in the southwest- 
northeast direction and has a counterpart at 1.3 cm, elongated 
also in the same direction. The spectral index, ~ 1, indicates 
that the emission is thermal and partially thick. Such a spec- 
tral index has been found for thermal radi oiets associated with 
sources powering molecular outflows (e . g.. iBeltran et al.ll2~001t 
Zapata et al. 2004; Gibb & Hoare 200 7]), and for hypercompact 
H 11 regions (e. g JCarral et all 1 9971; llgnace & ChurchweHl2004t 
lKurt3l2005h . There is an infrared source with strong infrared ex- 



cess associated with the centimeter source, suggesting that the 
massive YSO is likely surrounded by circumstellar material (see 
Fig. |2]i. In addition, the dust continuum emission at 1.2 mm 
consists of a main condensation whose peak is clearly associ- 
ated with the centimeter source (the offset is < 5", the posi- 
tional uncertainty of the IRAM 30 m) , as well as with H2O m aser 
and NH 3 emission dPalla et alJll99U iMolinari et al.lll996l) . The 
spatial coincidence of the millimeter peak with the centimeter 
source suggests that the underlying massive star is in a very 
young evolutionary stage, not having pushed away the surround- 
ing material yet. This is reinforced by the detection of a compact 
source at 7 mm at the position of the centimeter source. Note 
that 38% of the 7 mm flux density comes from the dust emis- 
sion from the envelope/disk associated with the massive YSO, 
and for this dust component we estimated a mass of ~ 13 Mq. 
Note also that the 7 mm emission, similar to the 3.6 and 1.3 cm 
emission, is slightly elongated in the so uthwest-nort h east d irec- 
tion. Regarding the outflow emission, IZhang et al.l (120051) de- 
tect no CO (2-1) outflow in the region using the NRAO 12 m 
Telescope. However, CO (1-0) emission had been previously de- 
tec ted using the IRAM 30 m with a linewidth of 6.1 km s _1 
bv lWouterloot & Brand ( 19891). with a sensitivity similar to that 



achieved by Zhang et al 



(2005). tions 
are required to confirm or discard the presence of an outflow as- 
sociated with the massive YSO, and thus to properly interpret the 
nature of the elongated emission at 7 mm. The values for the dust 
temperature and mass derived from the SED fit are around 30 K 
and 23 Mq, which are near the median values derived for the 
sources of this work. Finally, note that there are 2MASS sources 
showing infrared excess, surrounding the centimeter source and 
falling within the millimeter condensation (see Fig. [2), indicat- 
ing that a cluster of low-mass YSOs is possibly forming around 
the massive YSO. 



4.2.3. IRAS 18123-1203 

We detected four sources toward 118123 at 3.6 cm within the 
central 50", as well as a dust condensation consisting of differ- 
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Fig.l. IRAS 00117+6412. Top: Black contours: VLA 3.6 cm con- 
tinuum emission. Levels are -6, -3 and 3 to 33 in steps of 3, times 
0.038 mJy bearrT 1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -4, -3, 3, to 8 in steps of 1, times 0.067 mJy beam -1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 3 to 13 in steps 
of 1, times 8 mJy bearrT 1 . White contours: 3.6 cm emission as in top 
panel. Symbols in Figures [T] to [TO] are described here. Red ellipse indi- 
cates the error ellipse of the IRAS source. White and/or black crosses 
indicate the centimeter sources from Tables [6] and [7] White filled cir- 
cles indicate CH 3 OH masers, white filled triangles H 2 masers, and 
white filled squares OH masers. Color stars indicate 2MASS sources 
with the color and symbol corresponding to different infrared excess as 
described in § 13.31 Synthesized beams at the bottom-left corners cor- 
respond to the first wavelength indicated at the top of the panel, and 
beams at the bottom-right corners, when available, correspond to the 
second wavelength at the top of the panel. 



ent peaks (see Fig. [3j. However, neither the centimeter emis- 
sion nor the main millimeter emission are coincident with the 
IRAS source, and the centimeter and millimeter emissions are 
not overlapping. In addition, there are no 2MASS sources asso- 
ciated with the main dust condensation or the centimeter sources 
(except maybe for VLA 4). About 10" to the south of the 
main dust condensation there is a faint and compact millime- 
ter source at 3<x, associated with a 2MASS source with some 
infrared excess, and within the IRAS error ellipse. The faint 



Fig. 2. IRAS 04579+4703. Top: Black contours: VLA 3.6 cm con- 
tinuum emission. Levels are -4, -3 and 3 to 7 in steps of 1, times 
0.022 mJy beam -1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -4, -3, 3, to 5 in steps of 1, times 0.085 mJy beam -1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 12 in steps 
of 2, times 7 mJy beam -1 . White contours: VLA 7 mm continuum emis- 
sion. Levels are -4, -3, 3, to 7 in steps of 1, times 0.094 mJy beam - '. 
Color stars indicate 2MASS sources and white filled symbols refer to 
masers (see Fig.[T]for more details). 



JHK magnitudes of this 2MASS source, its marginal detection 
at 1 .2 mm, and the non-d etection of NH3 emission toward the re- 
gion (Molinari et al. 1996) suggest that the YSO associated with 
the IRAS source should be of low-mass. However, this is in con- 
flict with the assumed IRAS bolometric luminosity of 7900 Lq. 
Thus, either this object is peculiar, or its distance is overesti- 
mated. 

The spatial anticorrelation between the centimeter and the 
millimeter emission is somewhat puzzling. We first considered 
the possibility that the centimeter emis sion is tracin g externally 
ionized gas, and searched the catalog of Reed ( 20030 for nearby 
OB stars. We did not find any OB star to the northe ast of the cen- 
timet er emission which could account (following Lefloc het alj 
\199% for the total centimeter flux of VLA 2, VLA 4 to VLA 6 
(1.6 mJy), and thus we ruled out the possibility that these cen- 

3 Catalog available at http://othello.alma.edu/~reed/OBfiles.doc 
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Fig. 4. IRAS 18171-1548. Top: Black contours: VLA 3.6 cm con- 
tinuum emission. Levels are -5, -3 and 3 to 9 in steps of 2, times 
0.185 mjy beam -1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -4, -3, 3, to 6 in steps of 1, times 0. 1 12 mjy beam -1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 5 in steps 
of 1, times 16 mjy beam -1 . White contours: 3.6 cm emission as in top 
panel. Orange filled stars indicate IRAC-Spifzer sources classified as 
Class 0/1 sources (see § 13. 3t , and with no 2MASS counterpart. Color 
stars indicate 2MASS sources and white filled symbols refer to masers 
(see Fig.[T]for more details). 



Fig.3. IRAS 18123-1203. Top: VLA 3.6 cm map. The box shows 
the region zoomed in the other panels. Middle: Black contours: VLA 
3.6 cm continuum emission. Levels are -4, -3 and 3 to 7 in steps of 
1, times 0.043 mjy beam -1 . White contours: VLA 1.3 cm continuum 
emission. Levels are -4, -3, 3, and 4, times 0.063 mjy beam -1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 5 in steps 
of 1, times 15 mjy beam -1 . White contours: 3.6 cm emission as in mid- 
dle panel. Color stars indicate 2MASS sources and white filled symbols 
refer to masers (see Fig.Q~]for more details). 



timeter sources were externally ionized globules. Another possi- 
bility is that the centimeter emission traces extragalactic objects, 
suggested by the double radio source morphology. However, 
there are no sources detected at 21 cm in the NVSS (the spectral 
index between 21 and 3.6 cm being typical of thermal emission), 
and the spectral indices between 3.6 and 1.3 cm do not discard a 
free-free thermal nature for the centimeter emission (except for 
VLA 5). Thus, the centimeter emission to the north of the IRAS 
source could be produced by a group of low-mass YSOs driv- 
ing thermal radiojets, whose extended emission may have been 
resolved out at 1.3 cm. 
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Table 11. Parameters of the Spectral Energy Distribution fits, and dust thermal emission at 1.2 and 7 mm 









From SED Fit 


1.2 mm 


7 mm Dust Emission 
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free-free 
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S v h 




Region 


Source 


Spectral Index 
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(Mq) 




(mjy) 
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IRAS 00117+6412 


2-3 d 


+0.18 


1.5-1.7 


30-26 


21-30 


36 






IRAS 04579+4703 


1 


+ 1.1 


1.8-2.0 


32-28 


20-26 


23 


0.5+0.4 


8-19 


IRAS 18171-1548 


1 


-0.3 


1.7-2.0 


30-26 


100-170 


142 






IRAS 19045+0813 


7 


+0.4 


1.9-2.2 


29-26 


8-12 


10 






IRAS 22187+5559 


3 


+0.0 


1.9-2.1 


34-32 


10-20 


14 








5 


-0.06 


1.9-2.0 


32-31 


30-40 


32 






IRAS 22198+6336 


2 


+0.5 


1.7-1.9 


29-25 


20-30 


20 e 


1.2 + 0.7 


6-11 


NGC7538-IRS9 


2-3 d 


-0.10 f 


1.9-2.0 


35-32 


50-80 


66 s 







11 Dust and gas mass from IRAM 30 m 1.2 mm observations (TableO. 

b Flux density at 7 mm without the free-free contribution estimated from the 3.6 and 1.3 cm emissions. 

c Dust and gas mass from the 7 mm dust emission, assuming a dust temper ature in the 50-100 K range (Gomez e t al. 2003), a dust emissivity 
index of 2 and a dust mass opacity coefficient of 0.9 cm 2 g -1 at 1.2 mm l Ossenkopf & Hennindl 19941) . The higher value corresponds to a 
temperature of 50 K, and the lower to a temperature of 100 K. 

d Emission from the sources included in the range. 

c Estimated from the 850 yum flux density from Jenness et al] d 19951) . 

f Value adopted for the SED fit, due to the dete ction of the source at only one wavelength. 

s Estimated from the 1.2 mm flux density from Sandell & Sievers 120041). 



4.2.4. IRAS 18171-1548 

The centimeter source, detected at 3.6 and 1.3 cm, has a flat spec- 
tral index, a size of ~ 0.15 pc, and is clearly associated with an 
infrared source detected in the JHK bands from 2MASS and re- 
vealed by the IRAC-Spitzer images at 3.6-8.0 yum as a bright 
nebulosit y (but with no point source) emerging from a dark 
structure (Ku mar & Gravdl2007l) . Thus, the centimeter source is 
most likely tracing a compact H n region. Note that the compact 
H ii region lies in a zone with a low density of infrared sources 
(see the 2MASS /f-band image in Fig. 0), whose central part is 
coinciding with the dark structure revealed by IRAC. Note also 
that the IRAS source is slightly shifted to the west of the com- 
pact H ii region, and has a 2MASS source with strong infrared 
excess associated. The parameters derived for the compact H n 
region suggest that the underlying star is a Bl star, and from 
the fit to the SED we estimated a mass of the envelope around 
100 Mq for dust temperatures of ~ 28 K. The morphology re- 
vealed by the large-scale millimeter emission shows an elon- 
gated and clumpy condensation in the north-south direction co- 
inciding with the zone with a small number of 2MASS sources. 
Thus, the millimeter emission seems to trace the cloud of dense 
gas and dust in which the massive star was formed, which is the 
most massive cloud of this survey, as derived from the flux den- 
sity at 1 .2 mm. It is worth noting that the main dust condensation 
contains multiple peaks or subcondensations. Interestingly, the 
millimeter emission shows two peaks separated ~ 15", with the 
compact H n region lying in between them. This is suggestive of 
the compact H n region pushing out the surrounding material to- 
ward the north and south, with the cloud being in the process of 
disru ption, and decreasing its density. Thi s could explain why no 
NH 3 (lMolinarietalll996h . and CS (2-1) dBronfman et all 19961) 
has been detected toward the region. 

4.2.5. IRAS 19045+0813 

The centimeter emission toward 119045 presents two main com- 
ponents (see Fig. |5j, the UCHu region candidate, and a sec- 
ond component with extended emission to the northwest of the 



UCH ii region. The UCH n region, falling within the IRAS er- 
ror ellipse, is associated with a 2MASS source with infrared ex- 
cess and with an IRAC-Spitzer source dKumar& Gravel 120071) 
classified as Class 0/1 source from the ^T-[3.6] vs [3.6]-[4.5] 
an d l~3.6H4.51 vs [4.51-^ 5.81 diagrams (§|33V Ther e is CS (2- 
1) teronfman et aill 19961) and H2O maser emission (Pall a et all 
1 1 99 lb detected in the region. The 1 .2 mm continuum emission is 
associated with the UCH 11 region, with three IRAC sources clas- 
sified as Class 0/1 (§ 13.31 ), and with faint emission at 7 mm. The 
7 mm emission has a secondary peak toward the south, in the 
same direction as the extension appreciated in the 1.3 cm emis- 
sion. The 7 mm emission does not seem to have contribution 
from thermal dust emission; therefore, the centimeter emission, 
up to 7 mm, is most likely tracing an ionized wind, as suggested 
by the spectral index of 0.2-0.6. The dust temperature derived 
from the SED, 26-29 K, is similar to the median value of the 
survey, and the associated mass is the smallest, around 8-12 Mq. 



The extended centimeter emission to the northwest consists 
of a curved structure with different subcomponents, VLA 2 to 
VLA 6, and is not overlapping with the millimeter emission, but 
starts when the millimeter emission finishes. This is reminiscent 
of a bright-rimmed cloud, as the centimeter emissio n follows the 
border of the dust cloud. We searched the catalog of lReedl {2003) 
for nearby OB stars to the northwest of the millimeter emission, 
and the closest star was a B2.5 star at ~ 2.5° to the north-west 
(70 pc), which is to o far away to accou nt for the flux of VLA 2 
to VLA 6, following lLefloch et ail d!997l) . Given the double lobe 
morphology and the negative spectral indices of some of these 
sources (see Fig. [5] and Table [5), they could be tracing extra- 
galactic objects. However, there are no sources detected at 21 cm 
in the NVSS, and we estimated a lower limit for the spectral in- 
dex of VLA 2 to VLA 6 between 21 and 3.6 cm of > -0.03 
(see Table |9), typical of thermal emission. Thus, we favor the 
interpretation that these sources are thermal, and have not been 
detected at 1.3 cm because their emission is extended and thus 
(partially) resolved out by the interferometer. 
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Fig. 5. IRAS 19045+0813. Top: VLA 3.6 cm map. The box shows the 
region zoomed in the other panels. Middle: Black contours: VLA 3.6 cm 
continuum emission. Levels are -5, -3 and 3 to 1 1 in steps of 2, times 
0.054 mJy bearrr 1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -5, -3, 3 to 11 in steps of 2, times 0.065 mJy beam -1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 10, in steps 
of 2, times 9 mJy bearrr 1 . White contours: VLA 7 mm continuum emis- 
sion. Levels are -4, -3 and 3 to 5 in steps of 1 , times 0. 100 mJy beam -1 . 
Color (orange filled) stars indicate 2MASS (IRAC-Spitzer) sources and 
white filled symbols refer to masers (see Figs.[T]and|4]for details). 



Fig. 6. IRAS 20293+4007. Top: VLA 3.6 cm map. The box shows the 
region zoomed in the other panels. Middle: Black contours: VLA 3.6 cm 
continuum emission. Levels are -5, -3 and 3 to 15 in steps of 2, times 
0.058 mJy beam -1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -4, -3, and 3 to 7 in steps of 1, times 0.080 mJy beam -1 . 
Grey scale: K-band 2MASS image. Bottom: Grey scale and black con- 
tours: IRAM 30 m continuum emission at 1.2 mm. Levels are 3 to 6 in 
steps of 1, times 9 mJy beam -1 . White contours: 3.6 cm emission as 
in middle panel. Color stars indicate 2MASS sources and white filled 
symbols refer to masers (see Fig.Q]for more details). 
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4.2.6. IRAS 20293+4007 

At 3.6 cm, 120293 is dominated by one strong source (see Fig. [6]). 
This source has two subcomponents, VLA 1 (south) and VLA 2 
(north), and the centimeter peaks are not clearly associated with 
any 2MASS source, although a few fall within the extended cen- 
timeter emission. At 1.3 cm there is one strong source slightly 
resolved at the position of VLA 2. The spectral index suggests 
that VLA 2 is a thermal source, and is likely tracing the UCH n 
region associated with the IRAS source. However, the IRAS 
source could be also associated with a 2MASS source showing 
infrared excess that falls just at the eastern edge of the centimeter 
emission. The emission at 1 .2 mm is extended and clumpy, and 
shifted to the east with respect to the centimeter emission (sim- 
ilarly to 1001 17). There is no millimeter emission directly asso- 
ciated with the centimeter emission, suggesting that the UCH n 
region (possibly VLA 2) has disrupted and pushed the surround- 
ing gas away. This co uld be again the rea son why NH3 is not 
detected in the region (Molinar ietal.1119961) . 

4.2.7. IRAS 22187+5559 

122187 contains the most extended centimeter source of this sur- 
vey (see Fig. |7). The morphology of the centimeter emission 
reveals two main subcomponents, encompassed by a common 
structure extending toward the northwest. The strongest sub- 
component, VLA 5, located in the north, is compact, has an 
extension toward the northwest and could be associated with a 
2MASS source (with no infrared excess) lying at about 3" to 
the northwest from the centimeter peak. The second subcompo- 
nent, VLA 3, fainter and located in the south, is elongated in the 
southeast-northwest direction. At 1.3 cm, the emission strongly 
resembles the 3.6 cm emission, with three faint additional com- 
ponents within the extended emission, VLA 1 to the north of 
VLA 3, and VLA 2 and VLA 4 to the northwest of VLA 5. The 
spectral index of VLA 3 and VLA 5 is flat for both sources, 
while for the three faint sources, the spectral index varies from 
flat to +0.6. These three faint sourc es could be globule s exter- 
nally ionized by VLA 3 and VLA 5. lGarav et alld2006l) find in 
IRAS 16128-5109 a similar case of centimeter emission coming 
from externally ionized globules. 

The large-scale 1.2 mm emission also presents two main 
condensations, which are slightly (~ 7") displaced from VLA 5 
(millimeter peak to the east) and VLA 3 (millimeter peak to the 
south). Th e 1 .2 mm emission i s likely associated with the HCO + 
emission (Richards et al. 1987 ) and the CO emission detected 
in the region dWouterloot & Brandlll989l) . SCUBA maps at 850 
and 450 yum (|Pi Francesco et al. 120081) show condensations as- 
sociated with the 1.2 mm peaks (see Fig. O. We estimated a flux 
density of 0.7 + 0.2 Jy at 850 yum and 6.1 + 2.2 Jy at 450 yum for 
the condensation associated with VLA 3, and a flux density of 
1.9 + 0.2 Jy at 850 yum and 12.0 + 2.9 Jy at 450 /mi for the con- 
densation associated with VLA 5. In the mid-infrared there are 
two MSX sources whose morphology and position are very sim- 
ilar to the two 1.2 mm, 850 yum, and 450 yum condensations. In 
addition, the POSSII red plates reveal a zone of high extinction 
starting at the position of VLA 5 and extending to the east, which 
matches with the MSX and 1.2 mm sources (see Fig. [7]). Another 
zone of high extinction is also found to the south of VLA 3 (see 
Fig. |7J. These high extinction zones, traced by the millimeter 
and mid-infrared emissions, are likely tracing a dense molecu- 
lar cloud seen against a bright background. We note that while 
VLA 3 is associated with optical emission, VLA 5 lies within 
the molecular cloud (dark in the optical), and thus seems to be 



more embedded and probably younger than VLA 3. Note also 
that the extension of the centimeter emission toward the west of 
VLA 5 and VLA 3 is suggestive of the material surrounding the 
centimeter peaks having a lower density on the west side than 
on the east/southern sides, where the higher-density gas could 
be slowing down the expansion of the ionized gas. 

A cluster of 7 mm sources is associated with both VLA 5 
and VLA 3. The strongest 7 mm source, associated with VLA 5, 
is elongated in the east-west direction, surrounded by different 
fainter 7 mm sources, and possibly associated with the closest 
2MASS source, which is shifted only 2'.' 1 . For VLA 3, the emis- 
sion has several components with the strongest one located ~ 5" 
to the north of the centimeter peak. All the 7 mm emission as- 
sociated with VLA 3 and VLA 5 comes from thermal free-free 
emission. We built and fitted two different SEDs for VLA 3 and 
VLA 5, and obtained similar temperatures for both sources. The 
associated mass is higher for VLA 5 (~ 35 Mq) than for VLA 3 
(~15M Q ). 

4.2.8. IRAS 22198+6336 

This source i s loca ted in the L1204/S140 molecular complex. 
Tafal la et ail d 1993b . using the CS(l-O) transition, detect seven 
cores in a region of ~ 50', named A to G. L1204-G is asso- 
ciated with 122198. (Note that some authors name 122198 as 
L1204-A, which has the coordinates of IRAS 221 76+6303, so 
one m ust be cautious when using the labe ling of Tafalla et~al1 
Il993h . Observations of C O and C 18 Q by IJenness et alJ dl995h 
and observations of CO by IZhang et~ai] d2005l) show evidences 
of molecular out flow emission. There is also H2O maser emis- 
sion detec t ed bvlPalla et alJ dl99ll) . IValdettaro eTall d2002h and 
iFelli et ail d2007 ). who find several bursts of variable duration 
(200-500 days). 

Emission at 3.6 cm shows a compact single source lying 
within the IRAS error ellipse (see Fig. [HJ, with a counterpart 
at 1.3 cm, but with no 2MASS neither MSX sources associated. 
The centimeter peaks spatially coincide (~ 5" to the northeast) 
with a strong and comp act submillimeter co ndensation detected 
at 450 and 850 jt/m by IJenness et alJ d 1995b . In addition, emis- 
sion at 7 mm has been detected toward the source, with 44% 
of the flux at 7 mm coming from thermal dust emission. The 
mass of the 7 mm dust component is ~ 9 Mq. The source at 
1.3 cm and 7 mm is elongated in the southeast-northwest direc- 
tion (see Fig. [8]bottom panel), with a structure reminiscent of a 
radiojet, consistent with the spectral index measured in the cen- 
timeter range, ~0.5. Furthermore, the large-scale outflow found 
in this region is slightly elongated in the southeast-northwest di- 
rection. However, since a substantial part of the 7 mm flux den- 
sity is coming from thermal dust emission, at this wavelength 
we are probably observing the superposition of an ionized wind 
and a dust disk. If the direction of the outflow is the same as 
the direction of the elongation of the 1.3 cm and 7 mm emis- 
sions, one could infer from the measured deconvolved size of 
the 7 mm source, 3250 x 260 AU, an upper limit for the size of 
the disk, ~ 300 AU, as the disk is expected to be perpendicular 
to the 7 mm source. It would be necessary to study the outflow 
with high angular resolution to confirm its direction and thus 
the upper limit for the size of the disk. It is worth noting that 
as no near/mid-infrared sources are associated with the centime- 
ter emission (note that there is no distance effect because this is 
the source of the survey with the shortest distance), this source 
must be deeply embedded in cold dust (as indicated by the de- 
rived dust temperature of 25-29 K). We propose that VLA 2 is 
the youngest source of our survey and is tracing an early-type B 
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Fig. 7. IRAS 22187+5559. Top-left: Black contours: VLA 3.6 cm continuum emission. Levels are -5, -3 and 3 to 42 in steps of 3, times 
0.155 mjy beam -1 . White contours: VLA 7 mm continuum emission. Levels are -4, -3, and 3 to 9 in steps of 1, times 0.178 mjy beam -1 . Grey 
scale: K-band 2MASS image. Bottom-left: Grey scale and black contours: IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 12 in 
steps of 2, times 7 mjy beam -1 . White contours: VLA 1.3 cm continuum emission. Levels are -5, -3, 3, 5, 7 and 12 to 75 in steps of 7, times 
0.060 mjy beam -1 . Top-right: Grey scale: POSS-II red image. Black contours: 3.6 cm emission as in top-left panel. White contours: 1.2 mm 
emission as in bottom-left panel . Bottom-right: Grey scale and black contours: A-band MSX image (8.28 pm). White contours: SCUBA 850 pm 
emission. Levels are 3 to 17 in steps of 2, times 0.045 Jy beam -1 dDi Fr ancesco et al . 2008). Color stars indicate 2MASS sources and white filled 
symbols refer to masers (see Fig.fJJfor more details). 



protostar in an evolutionary stage similar to Class sources in 
the low-mass regime. 



4.2.9. NGC7538-IRS9 

We detected one source at 1.3 cm (at 3.6 cm we were limited 
by dynamic range), associated with the IRAS source and with a 
2MASS source at the northeastern border of an infrared nebu- 
losity (see Fig. [9]). The source at 1.3 cm splits up into two sub- 
components, VLA 2 and VLA 3, being VLA 3 clearly associated 
with the infrared source IRS 9. 

NH3 observations were carri ed out by Zheng et alJ d200ll) . 
and CO and H2 observations by iDavis et all dl998 | ). who dis- 
cover a H2 jet in the southeast-northwest direction. Iskndell et alJ 
(2005) detect three molecular outflows in the region: one HCO + 
high-velocity outflow associated with IRS 9 (VLA 3); another 
outflow associated with a protostellar core located ~ 25" to the 
southeast of IRS 9; and a third outflow, detected in the CO(l- 



0) transition, coincident with the H2 jet, w hose driving source i s 
traced by an infrared peak labeled "Al " bv lTamura et alJ d 199 lb . 
This infrared peak coincides with the position of the VLA 2 
peak. Thus, the centimeter source detected is probabl y the coun- 
terpar t of the infrared source driving the H2 jet dSandell et al.l 
2005). H2O, OH, and CH3OH maser e mission has been detected 
(iKameya et alJl99CilSandell et al.l2005l) in IRS 9 and in the out- 
flows. 

The large-scale millimeter emissi on toward NGC753 8 has 
been studied by a number of authors. lReid & Wil son ( 2005) map 
the submillimeter continuum emission at 450 and 850 pm, and 
detect a single source at both wavelengths, clearly associated 
with IRS 9, with a faint exte nsion to the south l i kely associated 
with the infrared nebulosity. San dell & Siev ers (2004) map the 
continuum emission at 350 pm and 1.3 mm, finding similar re- 
sults. 



Th e region was observed at 7 mm by Ivan der Tak & Mentenl 
d2005l) with the VLA in the C and A configurations. These au- 
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Fig. 8. IRAS 22198+6336. Top: Black contours: VLA 3.6 cm con- 
tinuum emission. Levels are -6, -3 and 3 to 21 in steps of 3, times 
0.027 mJy bearrT 1 . White contours: VLA 1.3 cm continuum emission. 
Levels are -5, -3, and 3 to 9 in steps of 2, times 0.063 mJy bearrT 1 . 
Grey scale: K-band 2MASS image. The box shows the region zoomed 
in the bottom panel. Bottom: Grey scale and black contours: VLA 7 mm 
continuum emission. Levels are -5, -3, and 3 to 7 in steps of 2, times 
0.226 mJy beairT 1 . White contours: 1.3 cm emission as in top panel. 
Color stars indicate 2MASS sources and white filled symbols refer to 
masers (see Fig.[T]for more details). 
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thors detect one source at the position of VLA 3 (positions agree 
within Of! 3), which is elongated in the north-south direction in 
the A-configuration observations. This is consistent with the de- 
convolved position angle that we have obtained toward VLA 3, 
10°. The present observations at 1.3 cm add a new point to 
the centimeter SEP , which was fitted through an ionized accre- 
tion flow model by van der Tak & Mentenl d2005l) . Our value at 
1.3 cm is about a fa ctor of 2 higher than that pre dicted by the 
model (see Fig. 4 of Ivan der Tak & Mentenll2005l) . most likely 
because our observations were carried out in the D configuration 
and thus picked up more extended emission. The dust temper- 
ature from the SED fitted by us (which includes both VLA 2 



Fig. 9. NGC 7538-IRS9. Top: Grey scale: 21 cm emission from NVSS. 
White contours: VLA 3.6 cm continuum emission. Levels are 3, 5, 10, 
15, 20, 50, 100, 200, 300 and 400, times 1.28 mJy beairT 1 . The box 
shows the region zoomed in the other panels. Middle: White contours: 
VLA 1.3 cm continuum emission. Levels are -5, -3 and 3 to 11 in 
steps of 2, times 0.204 mJy bearrT 1 . Grey scale: K-band 2MASS im- 
age. Blue lines correspond to molecular outflows ( Sandell et al. 2005). 
Bottom: Grey scale and black contours: IRAM 30 m continuum emis- 
sion at 1.3 mm. Levels are 5 to 35 in steps of 3, times 30 mJy bearrT 1 
lSandell&S 'ieversll2004h . White contours: 1.3 cm emission as in mid- 
dle panel. Color stars indicate 2MASS sources and white filled symbols 
refer to masers (see Fig.fT]for more details). 



16 



A. Sanchez-Monge et al.: Survey of cm and mm emission from massive star-forming regions 




Right Ascension (J2Q00) 

Fig. 10. IRAS 23448+6010. Top: Black contours: VLA 3.6 cm con- 
tinuum emission. Levels are -5, -3 and 3 to 13 in steps of 2, times 

0. 030 mjy beam 1 . White contours: VLA 1.3 cm continuum emis- 
sion. Levels are -5, -3, 3, 5, and 7, times 0.076 mJy bearrT 1 . Grey 
scale: K-band 2MASS image. Bottom: Grey scale and black contours: 
IRAM 30 m continuum emission at 1.2 mm. Levels are 2 to 5 in steps of 

1, times 1 1 mjy beam -1 . White contours: 3.6 cm emission as in bottom- 
left panel. Color stars indicate 2MASS sources and white filled symbols 
refer to masers (see Fig.Q]for more details). 



and VLA 3) using a modified blackbody law, is near the me- 
dian value of the survey, and the associated mass, 50-80 Mq, 
is among the highest in the survey. The small deconvolved 
size of the 1.3 cm source, its association with a 7 mm source 
( Ivan der Tak & Mentenll2005 | ) and with the peak of a compact 
1.3 mm source (ISandell & Sieversl l2004). the high infrared ex- 
cess of the 2MASS source associated with VLA 3, together with 
the fact that th e centimeter fluxes are consis tent with an ionized 
accretion flow (Ivan der Tak & Mentenl2 005) suggest that VLA 3 
(and probably VLA 2) is among the youngest sources in the sur- 
vey. 

4.2.10. IRAS 23448+6010 

The 3.6 cm emission toward this region shows a number of 
sources that lie in a region wi th a cluster of 2M ASS sources 
showing strong infrared excess (lKumaretalJl2006L see Fig. [TO]). 



The 1 .2 mm emission of the region is faint and extended and sur- 
rounds the northern and western edges of the centimeter emis- 
sion. Three of the centimeter sources, VLA 1 to VLA 3, have 
spectral indices compatible with thermal free-free optically thin 
emission, and thus could be UCH n regions, although with dif- 
ferent properties: while VLA 3 is the only centimeter source de- 
tected at 1.3 cm, and has no 2MASS counterpart, VLA 1 and 
VLA 2 are only detected at 3.6 cm and are clearly associated 
with 2MASS sources, with VLA 2 falling inside the IRAS error 
ellipse. As the 2MASS source coincident with VLA 2 is not de- 
tected in the /-band, we estimated its infrared excess from the 
( H — K) color alone, which i s 0.63, typical of T-Tauri stars (e. g., 
iKenyon & Hartmann II 19951) . and thus included in the group of 
sources with infrared excess < 0.4 (following § 13.31 ). We note 
that the centimeter emission of VLA 2 to VLA 6 is reminiscent 
of a ring-like structure and that this emission could be produced 
by the ionizing photons from the UCH n regions. However, none 
of these sources lie at the center of the ring, and in addition, the 
rate of ionizing photons coming from VLA 2 or VLA 3 is too low 
(more than one order of magnitude) to produce t he centimeter 
emiss ion of VLA 4, VLA 5 and VLA 6 (following Lefloc h et all 
Il997l) . Thus, the centimeter emission of 123448 may be tracing 
a young cluster of B-type stars, which could be dispersing the 
surrounding material and would be responsible for the morphol- 
ogy of the millimeter emission. iGarav et all ([2006) find a simi- 
lar case of a cluster of centimeter sources in IRAS 15502-5302, 
but with the millimeter emission still associated with the cen- 
timeter sources. HCO + and HCN em ission has been detected 
toward 123448 dRichards et all [l987l) as well as CO emission 
dWouterloot & Brandlll989l) . 

4.3. On the evolutionary stages of the sources of the survey 

In order to estimate the evolutionary stage of the sources pre- 
sented in this paper, we searched the literature for the presence 
of dense gas, outflow and maser emission. All this information 
is summarized in Table Q~2] The table shows that there are five 
sources with dense gas and maser emission associated, being 
thus possibly the youngest regions in the su rvey. This i s sup- 
porte d by the evolutionary classifications of iHill et all |2007) 
and lLongmore et all d2007l) for the formation of high-mass stars. 
According to these authors, the earliest stages corresponding to a 
high-mass "starless core" are characterized by cool (< 20 K) mil- 
limeter/dense gas emission with no maser or centimeter emission 
associated yet. After this stage, the first signposts of the mas- 
sive protostar are traced by warm (> 50 K) millimeter emissio n 
and class I/II-methanol maser emission (Elling sen et al J 120071) . 
which corresponds to the "hot-core" phase, not detected in the 
centimeter range yet. Finally, when the YSO starts the hydro- 
gen burning, it starts ionizing the surrounding gas and it can be 
detected in the centimeter range, while the YSO is still accret- 
ing part of its final mass. H2O and OH maser emission may be 
present at this last stage (Ellingsen et al. 2007). Thus, as a first 
approach, the UCH n regions associated with millimeter contin- 
uum emission, dense gas and maser emission can be considered 
younger than the UCH 11 regions without. 

In an attempt to better classify our UCH 11 regions in an 
evolutive sequence, we studied the morphology of the centime- 
ter and millimeter emissions, and the relation between them. 
From the separation between both peaks of emission, listed in 
Table [T2l one can see that there are four UCH 11 regions clearly 
associated with a compact millimeter peak (104579, 119045, 
122198, and NGC7538-IRS9), suggesting that the natal cloud 
has not yet been significantly disrupted by the massive YSO. 



A. Sanchez-Monge et al.: Survey of cm and mm emission from massive star-forming regions 



17 




10° 10 1 10 2 10 3 10 4 10 5 10 6 10° 10 1 10 2 10 3 10 4 10 5 10 6 



?i(|im) 



Fig. 11. Spectral Energy Distributions. Green dots: 2M ASS data. Yellow do ts: IRAC-Spitzer data. Blue dots: MSX data. Red d ots: IRAS 
data. Orange dots: millimeter and submillimeter data from lJenness et ID (QUI) for IRAS 22198+6336, and from Sandell & Sievers (2004) for 
NGC7538-IRS9. Black dots: this work, and submillimeter data from SCUBA (JCMT) for IRAS 22187+5559. Dashed lines: modified blackbody 
and free-free optically thin fits for the far-infrared to centimeter wavelengths. Black line: sum of the modified blackbody law and the free-free 
optically thin law. The results of the fit are shown in Table II II Note that there is no attempt to fit the near- and mid-infrared wavelengths. 
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Table 12. Evolutionary indicators of the sources, classified from smallest to largest disruption degree (see Section |4~3l 





VLA 


cm 


cm-mm 


mm 


Disruption 


Dense 




Maser Emission 




2MASS 


MSX 


Region 


Source 


Structure 


Offset 
(pc) 


Structure 


Degree 11 

(%) 


Gas 


Outflow 


CH3OH 


H2O 


OH 


IR Excess 


Emission 


IRAS 22198+6336 


3 


disk / radiojet 


£0.03 


compact 15 


7 b 


NH 3 , CS 


yes 




yes 


no 


_ 


no 


IRAS 19045+0813 


7 


compact 


£0.04 


compact 


46 


CS 




no 


yes 


no 


orange 


yes 


IRAS 04579+4703 


1 


compact 


£0.06 


compact 


56 


NH 3 


yes 1 ' 


no 


yes 


no 


red 


ves 


NGC 7538-IRS9 


3 


compact 


-0.09 


compact 13 


69 b 


NH 3 , HCN 


yes 


yes 


yes 


yes 


red 


yes 


tdac ioni t e AO 
LKAb 181/1 — 1548 


1 


compact 


~0.08 


peak+ext 


78 


no 




no 


no 


no 


yellow 


yes 


IRAS 22187+5559 


5 


peak+ext 


-0.10 


peak+ext 


78 


HCO+ 


yes d 




no 


no 


blue^ 


yes 


IRAS 00117+6412 


2-3 


compact 


-0.13 


peak+ext 


81 


NH 3 , CS 


no 1 


no 


yes 


no 


yellow 


yes 


IRAS 23448+6010 


2 


ext 


-0.19 


dispersed 


82 


HCN, HCO + 






no 


no 




yes 


IRAS 22187+5559 


3 


peak+ext 


-0.10 


peak+ext 


87 


HCO+ 


yes d 




no 


no 


green e 


yes 


IRAS 20293+4007 


2 


compact 


-0.25 


dispersed 


91 


no 






no 


no 


green e 


yes 



a Defined in Equation^ Th e error of the d isruption degree is - 4%. 

b 122198 850 }im map from |Jenness"eTaT] fl995]) , and NGC 7538-IRS9 1.3 mm map from Sandell & Sievers 1 2004). 
c Source not detected in the 2MASS catalog. 
d Weak outflow emission. 

c 2MASS sources falling inside the centimeter emission, but slightly offset from the centimeter peak. 

1 Outflow emission detected toward a source embedded in the dust condensation but not toward the centimeter source. 

B See Section flXiol 

Note.- The values of this table come from this work and from different references for each source. References for 1001 17: 4,6, 14, 17, 28,29; 104579: 7, 17, 26, 27, 28, 29; 118171: 4, 
17, 19, 25; 119045: 3, 4; 120293: 17, 19; 122187: 7, 9, 19, 20, 28; 122198: 7, 8, 10, 12, 15, 17, 19, 23, 24, 26, 28, 29, 31; NGC7538-9: 2, 11, 13, 16, 21, 30; 123448: 5, 7, 18, 19, 20, 
27,28. _ ^ ^ ^ ^ ^ ^ 

References: (l)IBaclnller et all 1991 (2 ) IBoogert et aU2002t (31IBrand et alfT99i (4)IBronfman et aU1996L (5)ICasoli et aU19861 (6)ICesaroni et all 19881 (7)IEdris et aU2007t (8) 
Furuya et al. 2003, (9) Gait et al. 1989, (10) de Gregorio-Monsalvo et al. 2006, (1 1) Hasegawa & Mitchell 1995, (12) Jenness et al. 1995, (13) Kameya et al. 1990, (14) Kim & Kurtz 
2006, (15) Larionov et al. 1999. (16) Mitchell & Hasegawa 1991, (17) Molinari et al. 1996, (18) Nguyen-Q-Rieu et al. 1987, (19) Palla et al. 1991, (20) Richards et al. 1987, (21) 
Sandell et al. 2005, (22) Slysh et al. 1994, (23) Tafalla et al. 1993, (24) Valdettaro et al. 2002, (25) van der Walt et al. 1995, (26) Wilking et al. 1989, (27) Wouterloot et al. 1993, (28) 
Wouterloot & Brand 1989, (29) Zhang et al. 2005, (30) Zheng et al. 2001, (31) Felli et al. 2007. 



We quantified the degree of disruption of the natal cloud by 
comparing the millimeter flux density encompassed in a region 
very close to the centimeter peak position with the total millime- 
ter flux density. More precisely, on one hand, the total millimeter 
flux density was evaluated as the flux density inside a 0.5 pc- 
diameter region (the typical size of massive d ense cores, e.g., 
iBeltran et alj|2006bl: IZinnecker & Yorkdl2007l) . centered on the 
centimeter peak position, S 0.5. The fixed diameter of 0.5 pc al- 
lows us to discard the contamination of the measured flux den- 
sity from nearby YSO or dense clouds. On the other hand, we 
evaluated the millimeter flux density inside a 0.2 pc-diameter 
region (the size of the IRAM 30 m beam at the distance of the 
farthest source of this survey) centered on the position of the 
centimeter peak, S 0.2. We defined the disruption degree as 



Disruption degree = 1 - . 

S 0.5 



(1) 



Thus, a low disruption degree indicates that the millimeter emis- 
sion is concentrated inside a 0.2 pc-diameter region around the 
centimeter peak position, while a high disruption degree in- 
dicates that the millimeter emission spreads out in a 0.5 pc- 
diameter region. Taking into account the pointing errors of the 
IRAM 30 m telescope, we estimated an error of ~ 4% for the 
values of the disruption degree. 

The disruption degree defined in this way should evaluate 
to what extent the massive YSO has pushed out and dispersed 
the natal gas and dust. In Table 1 121 the regions are ordered from 
smallest to largest disruption degree, listed in column (6). From 
the table, one can see that the sources with smallest disruption 
degrees, < 70%, are 122198, 119045, 104579, and NGC 7538- 
IRS9. If the cloud is disrupted by the expansion of the UCH 11 



region then one would expect to find a correlation between the 
size of the centimeter source and the disruption degree. Such 
a correlation is clear in Fig. Q~2j where we plot the disruption 
degree versus the size of the centimeter source. 

It should be mentioned that the disruption degree and the 
centimeter size could be affected not only by the expansion of 
the Hn region (an evolutionary effect ), but also by the i nitial 
density of the natal cloud. For example, de Pree et all ( 1 19951) find 
that for initial higher densities of the natal cloud the sizes of the 
UCH 11 regions associated are smaller. Similarly, the disruption 
degree could be affected by the initial density of the cloud: the 
smaller the initial density, the larger the disruption degree. In 
addition, inhomogenities in the medium surrounding the mas- 
sive YSO could produce an accumulation of matter where the 
density is higher, shifting the peak of millimeter emission from 
the centimeter peak, as could be the case of 122187 (see § 14. 2. 7b . 
However, several arguments are against the interpretation of the 
disruption degree as tracing only the initial conditions. First, in 
the first stages of the formation of a massive YSO (assuming 
no close companions), the YSO has to be associated with mil- 
limet er flux at least as strong as its surroundings (e. g. jHill et alj 
120071) . Second, we found small values of the disruption degree 
for the regions which we previously classified as less evolved 
judging from dense gas and maser emission (see above). In ad- 
dition, we note that our sample does not show an important cor- 
relation between the size of the centimeter source and the bolo- 
metric luminosity. So, the disruption degree seems to be tracing 
more the evolution of the massive YSO than the initial density 
of the cloud. 

Even more relevant is that our classification from the disrup- 
tion degree is consistent with the infrared properties and in par- 
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Fig. 12. Disruption degree of the natal cloud defined as one minus the ratio of the millimeter flux density picked up within a ~ 0.2 pc-diameter 
region (the size of the IRAM 30 m beam at the distance of the farthest source of this survey) centered on the centimeter peak position and the 
millimeter flux density picked up within a ~ 0.5 pc-diameter region (the typical size of massive dense cores) also centered on the centimeter peak 
position, versus the size of the centimeter source (from TablellOt. For the disruption degree the errors have been estimated taking into account the 
pointing errors of the IRAM 30 m telescope. For the cm-size we estimated the error from the gaussian fit obtained with the task JMFIT in AIPS. 



ticular with the infrared excess derived from the 2MASS data, 
which should also be indicative of the evolutionary stage (§ 13. 3b . 
We found that the source with the smallest disruption degree in 
our survey, 122198, has no near- (2MASS source) neither mid- 
infrared (MSX source) emission, suggesting that this YSO is 
deeply embedded in gas and dust, which is opaque to the in- 
frared emission. Regions with (slightly) larger disruption degree 
are 119045, 104579, and NGC 7538-IRS9, the three showing red 
or orange infrared excess (§[33}. For these sources the centime- 
ter emission is compact and spatially coincident with the mil- 
limeter emission, which is compact as well. Sources with larger 
disruption degrees as 118171 and 1001 17, for which the millime- 
ter emission appears clumpy, extended, and/or displaced from 
the centimeter peak, have smaller infrared excesses (yellow). 
Finally, sources with the largest disruption degrees (123448, 
120293) have the centimeter emission less compact, with the mil- 
limeter emission dispersed around the centimeter source, and 
the infrared excess being among the smallest in the sample. 
Thus, our sample contains intermediate/high-mass star-forming 
regions which can be ordered in an evolutionary sequence, from 
regions containing massive YSOs deeply embedded in the natal 
cloud, to regions with massive YSOs whose ionizing photons are 
in the process of completely disrupting the natal cloud. 

5. Conclusions 

We observed, with the IRAM 30 m telescope and the VLA, 
the continuum emission at 1.2 mm, 1.3 cm and 3.6 cm to- 
ward a sample of 11 intermediate/high-mass star-forming re- 
gions not previously observed at millimeter and/or centimeter 
wavelengths. Nine were observed in the millimeter range and 10 
in the centimeter range (8 coinciding with the regions observed 
at IRAM 30 m). Additionally, four of the most promising regions 
were observed with the VLA at 7 mm. Our main conclusions can 
be summarized as follows: 



1 . From the IRAM 30 m observations we found 1 .2 mm emis- 
sion in all regions, but 118212. Four regions have very ex- 
tended and weak emission, while the other six show a clear 
emission peak with some substructure. We derived the mass 
of gas and dust traced at 1.2 mm, finding values from 10 to 
140 Mq. 

2. From the centimeter VLA observations we found 3.6 cm 
emission in all regions, being most of them from compact 
sources. At 1.3 cm, all regions, but 118123, have emis- 
sion associated with some 3.6 cm sources and a centimeter 
source is found within or near the IRAS error ellipse. We 
estimated the spectral indices of the central sources (those 
lying within the central 50" around the IRAS source) be- 
tween 3.6 and 1.3 cm. Additionaly, we searched the NVSS 
for emission at 21 cm, and estimated the spectral index be- 
tween 21 and 3.6 cm. For most of the sources associated 
with the IRAS source we found flat or positive spectral in- 
dices, suggesting that the emission comes from compact or 
UC H ii regions, thermal radiojets or externally ionized glob- 
ules. Deconvolved sizes of the centimeter sources range from 
<0.01 to 0.3 pc, and the physical parameters of the Hn re- 
gions indicate that the ionizing stars are early B-type stars. 

3. For the four regions observed at 7 mm the emission was par- 
tially resolved. Taking into account the spectral index be- 
tween 3.6 and 1.3 cm, we were able to derive the dust con- 
tribution to the 7 mm continuum emission for 104579 and 
122198, which is of about 38% and 44%, respectively. 

4. By combining our data with the infrared surveys of 2MASS, 
MSX, and IRAS, and with IRAC-Spitzer data (when avail- 
able), we built the SED and fitted a modified blackbody law 
taking into account the contribution of a free-free optically 
thin law for the centimeter observations. We found median 
values of ~ 25 Mq for the mass, ~ 29 K for the dust tempera- 
ture, and ~ 1 .9 for the dust emissivity index for the envelope 
of each intermediate/high-mass YSO. 
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5. By studying the morphology of the centimeter and millime- 
ter emission, and the separation between both peaks of emis- 
sion, we found a correlation between the degree of disruption 
of the natal cloud (estimated from the millimeter flux picked 
up within a ~ 0.2 pc-diameter region, centered at the posi- 
tion of the centimeter peak, and the millimeter flux picked up 
within a ~ 0.5 pc-diameter region), and the size of the cen- 
timeter source. This allowed us to establish an evolutionary 
sequence for the UCH n regions studied in this work, which 
is consistent with the expected evolutionary stage from the 
dense gas, outflow, and maser emission, as well as with the 
infrared properties of the sources associated with the UCH n 
regions. Thus, we classified our sources from the youngest, 
122198, similar to the Class low-mass sources, with no 
near- neither mid-infrared emission, to the most evolved, 
123448 and 120293, with the centimeter source(s) having dis- 
rupted the molecular natal cloud. 
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